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ABSTRACT 

The standard model for extragalactic variable radio sources comprises an 
Isotropically expanding plasmoid with "frozen" magnetic flux and an electron 
distribution which evolves adlabatlcal ly. This model leads to the following 
relation between the peak luminosity m and the relevant frequency v m , , 
which are functions of time: L « v N where N a (7n + 5)/(4n + 5). In 

v ,m m 

this expression, n is the spectral index in the optically thin part of the 
spectrum, where L y * v“ n . For n in the range 0.5 to 1.5, the standard model 
yields N in the range 1.2 to 1.4. By contrast, analysis of observational data 
yields estimates of N in a small range about the mean value 0.4, in clear 
contradiction with the standard model. 

The model is here modified to comprise finite flux tubes, either rooted 
in a parent object or forming closed toroids. The expansion rate along the 
flux tube may differ from the expansion rate transverse to the flux tube. The 
electron distribution is assumed to remain isotropic and to evolve 
adiabatically. This model yields a good match to the above observational 
relation in the case that the flux tubes expand along their length but not in 

t 

the transverse direction. 

Subject Headings : radio sources: galaxies - quasars - radio emission 

^Also Department of Applied Physics, Stanford University 
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EXTRA6ALACTIC VARIABLE RADIO SOURCES 
I. Introduction 

It has been known, since the observations of Dent (1965), that radio 
emission In the millimeter and centimeter range from quasars and some radio 
galaxies Is variable on a time scale of years, months, and sometimes shorter 
Intervals. Soon after these variations were discovered, it was suggestd by 
Shklovsky (1965) that the varying components are sources of synchrotron 

radiation which are initially optically thick but, as they expand, become 

optically thin at progressively longer wave lengths. 

This proposal was analyzed in more detail by van der Laan (1966), who 

analyzed the properties of a model comprising a uniform and spherical cloud of 
magnetic field and relativistic electrons. The electron distribution was 
assumed to be Isotropic, with a power-law dependence on energy. The magnetic 
field was assumed to be "frozen" in the cloud so that magnetic flux is 
conserved. The expansion was assumed to be spherically symmetric so that 
there is no change of topology of the magnetic field and the electron 
distribution remains Isotropic. A further assumption was that the 
relativistic gas cools adiabatical ly during expansion. 

This model was compared with observational data by Kellermann and 

Pauliny-Toth in 1968, and the agreement was quite satisfactory. However, the 
situation has changed in recent years. 

Andrew, MacLeod, Harvey and Medd (1978) have analyzed results of a ten 
year study of extragalactic variable sources at centimeter wavelengths. They 
examined 99 variable sources and found that bursts of particle acceleration 
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occur In the average variable source about once every one and one-half to two 
years. 

At any time In the evolution of the source, the flux (per unit frequency) 
F has a maximum value F %1 _(t) at a frequency v_(t). Andrew et al. 

V v fill mi 

(1978) found that, for any source, these quantities vary In such a way that 
they satisfy approximately the relation 

'W 1 ’ - • »•») 

They found that observational ly determined values of N fall in a narrow range 
centered on 0.4: 


N 


0.4 


+0.15 

- 0.2 


( 1 . 2 ) 


By performing a similar analysis of data compiled by Altschuler and Wardle 
(1976, 1977), Andrew et al . (1978) determined N values in the range N = 0.4 + 
0.15. 

According to the SVDL (Shklovsky-van-der-Laan model), N should be related 
to n, the power-law index of the radio emission, by the relation 

7n + 5 

N = . (1.3) 

4n + 5 


Hence, for n = 0.5, 0.75, 1.0, 1.5, we expect that. N = 1.21, 1.28, 1.33, 1.41, 
respectively. It is clear that this prc":rty of the SVDL model is 


incompatible with the observational data summarized by Andrew et al . (1978) 
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II. Flux-tube Model 

Since the properties of the SVDL model are broadly consistent with 
observational data. It seems reasonable to explore minor modifications of this 
basic model. The simplest variant which came to mind Is that the expansion Is 
anisotropic with respect to the magnetic field. This is the possibility which 
will be explored in this article. 

If one Is to consider the possibility that the expansion is anisotropic 
with respect to the magnetic field, it Is clear that one should also consider 
the possibility that the electron distribution is anisotropic. Discussion of 
this possibility is deferred for a subsequent article. For the time being, we 
consider that the electron distribution remains isotropic, due either to 
particle-particle scattering or to scattering by magnetic fluctuations, or to 
some other process. 

We consider a single magnetic flux tube of uniform radius R and length 
L. The flux tube may be rooted in an accretion disk, in which case it is to 
be assumed that the detailed geometry is such as to provide mirror action at 
the ends of the tube (Fig. 1(a)). Alternatively, the flux tube may be in the 
form of a toroid, in which case L measures the perimeter (Fig. 1(b)). As an 
additional simplification, we assume that the transverse component of the 
magnetic field (transverse to the major direction measured by L) is small 
compared with the longitudinal component of the magnetic field (in the 
direction measured by L). 

We assume that L and R vary with time t according to relations 

= Ljf*, R = R x f 


L 


( 2 . 1 ) 
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where f(t) is an arbitrary function of time and , Rj, etc., are the values 
of L, R, etc., at t = tj, assuming that f(tj) * 1. For unaccelerated 
expansion, f(t) = t/t ^ • For any reasonable assumption concerning the 

conductivity of the plasma, we find that the total flux of the tube Is 
constant so that the magnetic field strength B varies as follows, 

B = B^" 2 * . (2.2) 


For a relativistic gas with an isotropic pressure p, the pressure varies 
with density n according to the relation 


p « 



(2.3) 


The pressure of electrons with energy In a small range about the value E (eV), 
is proportional to nE. Hence, noting that 


n 



-U+2p ) 


(2.4) 


we see that 

E „ f "( 1/3) (X+2p ) 


(2.5) 


We suppose that the total number of relativistic electrons In the plasmoid has 
a power-law distribution of the form 


dN 


KE" m dE . 


( 2 . 6 ) 
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Since the total number of electrons Is fixed, and since the evolution of 
energy with time Is taken to be the adiabatic form (2.5), we find that 


K(t) 


f - ( 1 / 3 ) (a+2p ) (m-1) 


(2.7) 


If the lower and upper limits of w .he power-law spectrum (2.6) are taken to be 
E[_(t) and Ey(t), respectively, we see from Eq. (2.5) that 

r _ r f“( 1/3) (X + 2p ) p _ p f -(l/3)(X+2p) 

E L - E LJ f , E 0 - E ua f . (2.8) 


One may verify from (2.6), (2.7) and (2.8) that the total number of electrons 
Is fixed, being given by 


N 


T 



rr-(m-l) 

1 L.l 


p-(m-l) 1 

L U,1 J 



r-(m-l) 

b L.l 


(2.9) 


III. Radiation 

We now consider the synchrotron radiation produced by the electron 
distribution proposed In § II. For simplicity, we represent the radiation 
spectrum of a single electron S^(erg s^Hz"*) by the delta-function 
approximation 


S. = 2A.E 2 B 2 5(v-v ), A. - 10" 26,1 
1 ,v 1 p i 


(3.1) 


where 
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v p - A 2 E 2 B, A 2 - 10“ 5,3 . (3.2) 

On combining these equations with Eq. (2.6), we obtain the standard result 
that the total luminosity spectrum L v (erg s~*Hz’*), defined by 

L v * / dN S l,v * (3 - 3 > 

Is given by 

L v = A i A 2’ lKBn+lv " P • (3.4) 

where 

n = i (m - 1) . (3.5) 

The lower and upper bounds of this power-law spectrum are given by 

v L - A 2 E 3 B, V(j =■ AgE^B . (3.6) 

Following van der Laan (1966) and Scheuer (1967), we estimate the 

radiation In the self-absorbed part of the spectrum by noting that the 

brightness temperature cannot exceed the equivalent temperature of electrons 

chiefly responsible for the synchrotron radiation at any given frequency. On 

using the Raylei gh-Jeans formula for the self-absorbed form of the surface 

-2 -1 1 

emisslvlty (erg cm fc s A Hz -1 ) , 
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,SA ■ ¥ v • *3 " 10 


-32.0 


(3.7) 


and using Eq. (3.2) to express E In terms of v, we obtain the relation 


F v,SA ' V 1 / 2 v 6/!l . A 4 - Uf 30 - 3 . 


(3.8) 


If L and R remain comparable In magnitude for all time, we are 
effectively considering Isotropic expansion, that Is, we are reconsidering the 
SVDL model. It Is not possible that L should be much less than R. Hence, we 
consider only the case that L » R. With this Inequality, the effective area 
of the plasmold is effectively 2*RL, so that the self-absorbed part of the 
luminosity spectrum Is given by 

L v.SA * A 5 LRB- 1/2 v 5/2 , A. - lO' 29 - 5 . (3.9) 

We assume that the actual spectrum is given by (3.9) up to the frequency 

for which expressions (3.4) and (3.9) are equal, and by expression (3.4) 

for values of v above v m . We then find that the frequency v m of peak 
luminosity Is given by 

1 

v = [ A. A5 _1 A" 1 K I.'- R" 1 B n+3/2 l 0+5/2 , (3.10) 

m L l c a J 


v ,m 


and the peak luminosity L 


Is given by 
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v ,m 


r A 5/2 A (5/2)(n-l) A n K 5/? L n R n B 2n+5/2 1 7T?572 . 
1 1 2 5 i 


(3.11) 


On using Eqs. (2.1), (2.2), (2.7) and (3.5), we find that v (t) and 

m 

L v m (t) depend on time, through the function f(t), as follows: 


v 


m 



-(2/3)n+l . 

' n+5/2 x 


(10/3)n+4 
" nT5>2 p 


(3.12) 


and 




(19/3)n+5 
" n+5/2 p 


(3.13) 


where 


m,l 


[AjAj^A-ik^Lj'B ^ 3 / 2 ]"* 572 


(3.14) 


and 


1 

. r A 5/2 A (5/2)(n-l) A n K 5/2 L n R n R 2n+5/2 1 n+5/2 
v,m,l L A ! a 2 *5*2 L iri J 


(3.15) 


We see from the above equations that there should be a power-law relation 

between L and v , which we write as 
v ,m m 

• ( 3 - 16 ) 


where 
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2nx + (19n ♦ 15)p 

N - It _ ■ (3.17) 

A,p (2n + 3)x + (lOn + 12)p 

IV. Discussion 

We may readily verify that the case x • 1, p ■ 1, corresponding to 
Isotropic expansion, yields 



7n + 5 
4n + 5 


(4.1) 


This reproduces the familiar relationship between L , and v for the SVDL 

v ,m m 

model (van der Laan, 1966). 

We now consider two other special cases. The case \ ■ 0, p ■ 1 corre- 
sponds to a toroidal plasma for which the perimeter Is fixed but the minor 
radius grows progressively. It is clear that such a model would eventually 
break down, since the radius would at some stage exceed the perimeter. 
Nevertheless, we note that this model yields the following value of N: 



19n + 15 
lOn + 12 


(4.2) 


We also consider the case that X * 1, p ■ 0, corresponding to toroidal 
plasma of constant cross section, which Increases progressively In length 

2n 

Vo 


2n + 3 


(4.3) 
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We list the values of the above formulas for a range of values of n In 
Table 1. It Is clear that the second c*se, N 0(1 , provides an even worse fit 
to observational data than does the SVDL model, N^. On the other hand, the 
third model, characterized by N^g, appears to provide a quite acceptable fit 
to observational data. We see that for N In the range 0.65 < n < 1.50, N Is 
In the range 0.3 < N < 0.5. 

It Is clearly of Interest to pursue the case that X ■ 1, p ■ 0. For 

this case, Eq^. (3.12) and (3.13) yield the following dependence of v and 

m 

L m on f(t): 
v ,m 


v 


m 



(4.4) 


and 


where 


and 


v ,m 


L v ,m,l f 


-A 


4n + 6 


6n + 15 


4n 


6n + 15 


(4.5) 


(4.6) 


(4.7) 


Values of these Indices are given, for a range of values of n, in Table 2. We 
see that A varies only over the range 0.44 - 0.52 as n covers the range 0.5 
- 2.0, whereas A varies from 0.11 - 0.30 over the same range of values of n. 
Clearly, one cannot Infer from these relations the expected dependence of 
and l_ v ^ on t unless one can determine the expected dependence of f upon 
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t. If It cm be argued that L Increases linearly with time, so that 

f(l) « t then we see from Table 2 that, to fair approximation, v * t 1 ^ 2 . 

m 

As tests of the proposed model, It Is desirable to examine cases for 
which Individual bursts are well Identified to try to check the dependence of 
N upon n, given by Table 1, and the dependence of v |n and ^ on t, cjlven by 
Table 2. 

Finally, we Inquire Into mechanisms wi ,.i might produce plasmolds which 
expend progressively In length, while the cross section remains substantially 
constant. We first consider the property that the radius of the flux tubes 
does not change, or changes only slightly, during the evolution of the 
plasmold. This Is equivalent to the requirement that the mean magnetic field 
strength remains constant. If the dominant stress In the plasmold Is the 

magnetic stress, this condition Is equivalent to constant pressure, Implying 

that the plasmold Is moving In a constant-pressure medium. 

There are at least three processes which might lead to a progressive 
Increase In the length of the flux tube, as Indicated schematically In Fig. 
2. If the flux tube Is rooted In an accretion disk, then differential 
rotation of the disk may lead to a progressive stretching of the loop, as 

Indicated In Fig. 2(a). If, on the other hand, the plasmold comprises a 
toroid ejected from the disk, the steady Increase In length may be due to an 
Initial velocity of expansion of the toroid (Fig. 2(b)). This expansion could 
possibly be due to centrifugal motion, If the formation of the toroid Is 

accomplished by flare action and If, as Is the case In solar flares, the flare 
process leads to the formation of a dense flare plasma In the plasmold created 
by the flare. However, the centrifugal force would decrease with time more 
rapidly than the radially Inward force due to magnetic tension, so that the 
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expansion would ultimately be reversed. If the toroidal plasma is created 
with a high velocity in the radial (major radius) direction, we again find 
that magnetic tension would eventually reverse the expansion. 

A third possibility which comes to mind is that the plasmoid is 
turbulent, and that - as is well known - the turbulence leads to progressive 
stretching of magnetic field lines. If the external pressure is substantially 
constant, the magnetic field strength will remain constant. For any portion 
of the total flux tube comprising the plasmoid, the behavior is as described 
in this article. However, the total geometry is quite different: the overall 
shape may be spherical and the overall expansion may be spherically 
symmetric. Hence the model would be similar to the SVDL model except that the 
magnetic field strength remains constant rather than decreasing with time. 
This model is currently being examined. 
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Table 1 

Dependence of N upon n 
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Table 2 



Dependence of 

r and a on n 

n 

0.5 

0.44 

0.11 

0.75 

0.46 

0.15 

1 

0.48 

0.19 

1.5 

0.5 

0.25 


2 


0.52 


0.3 
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Fig. 
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Figure Captions 

. (a) Schematic model of flux tube rooted in accretion disk. 

(b) Schematic toroidal flux tube ejected from accretion disk. 

(a) Differential rotation in accretion disk leads to stretching of 
flux tube. (b) Schematic representation of toroidal flux tube, 
ejected from disk, that expands in major radius but not in minor 


radius 
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